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Abstract. We consider the so-called G-equation, a level set Hamilton-Jacobi 
equation, used as a sharp interface model for flame propagation, perturbed by 
an oscillatory advection in a spatio-temporal periodic environment. Assuming 
that the advection has suitably small spatial divergence, we prove that, as the 
size of the oscillations diminishes, the solutions homogenize (average out) and 
converge to the solution of an effective anisotropic first-order (spatio-temporal 
homogeneous) level set equation. Moreover we obtain a rate of convergence 
and show that, under certain conditions, the averaging enhances the velocity 
of the underlying front. We also prove that, at scale one, the level sets of the 
solutions of the oscillatory problem converge, at long times, to the Wulff shape 
associated with the effective Hamiltonian. Finally we also consider advection 
depending on position at the integral scale. 



1. Introduction 

We study the limit, as e — > 0, of the solution to the level-set equation 

r (i) nf = lZ)^x-I + (y(f,i),Dn-) in x (0,T) 

\ {ii) = no on X {0} . ^ ' ' 

Equation (|l.ip (i) is referred to as the G-equation, and is used as a model for 
flame propagation in turbulent fluids ( |23l |2l] ) • In that setting, the level sets of 
the function u"^ represent the evolving flame surface and —V is the fluid velocity 
field. At points where u'^ is differentiable and \Du''\ ^ 0, the level sets of u"^ move 
with normal velocity 

z.= l-(y(^,^),n) , 

where h = —Du'' /\Du''\ is the exterior normal vector of the front. When 1/ = 0, 
level sets move with constant speed sl = 1, which is called the laminar speed of 
flame propagation. 

We assume that the vector field V € C°'^(M^+^; M^) is Z^+^-periodic in 
both X and t, i.e., for all {x,t) G M^+i, A; G and s G Z, 

V{x + k,t + s) = Vix,t). (1.2) 

Our first result says that there exists a positively homogeneous of degree one, 
convex and continuous Hamiltonian H such that, as e — >• 0, the u^'s converge 
locally uniformly in X [0, oo) to the solution u of the initial value problem 
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ut = H{Du) in R^x[0,oo), 

(1.3) 

u = uq on X {0} . 

Although V is bounded, we do not assume that \V\ < 1, and, hence, the 
Hamiltonian 

H{x,t,p) = \p\ + {V{x,t),p) 

is not coercive in \p\ at every point {x,t). This lack of coercivity is the main 
mathematical challenge in the analysis. If either \V\ < 1 or the nonlinearity \Du\ 
were replaced with |-Dnp, then H would be coercive in \p\ and the problem would 
be within the scope of the theory developed in |19j . There are, however, relatively 
few homogenization results about noncoercive Hamiltonians [U El El Ell EZ] , and 
none of them deals with the particular structure considered here. 

The following simple example shows that, in the absence of coercivity, some 
additional assumption about the divergence of V is necessary in order for the n'^'s 
to have a local uniform limit. To this end, let V = V{x) be a smooth Z^-periodic 
vector field such that, in the cube Qi = (— ^, ^)^, ^(^) = — lOx, if |x| < 1/6, and 
V{x) = if |x| > 1/3. It is known that the u^'s have the control representation 

u'{x,t) = snp{uo{X.^{t))), (1.4) 

where the supremum is over ah functions G W^'°^{[0, t];R^) such that X^{0) = 
X and X'^(s) = k{s) + V{Xx{s)) with the controls k(-) satisfying \k\ < 1. If 
uq{x) = (p,x) with \p\ > 0, we see easily that lims^o\u'^(0,t)\ = for all 
t > 0. However, liminf^^Q u^(xe,t) > if t > and {xe}e is any point satis- 
fying |xe| = e/2. Roughly speaking, the problem with such a vector field V is 
that it traps the trajectories which start at the lattice points. If the divergence 
of V is sufficiently small, however, it is reasonable to expect that the controls are 
strong enough to overcome such traps. 

We assume that V has "small divergence", in the sense that, for all t € M, 

1 

a{t) = lldiv^y(-,f)||iiV(Q ^ >0 and a* = I a(s)ds>0, (1.5) 

ci Jo 

where c/ is the isoperimetric constant in the cube Qi (see, for instance, [E]), i.e., 
the smallest constant such that, for all measurable subsets E of Qi, 

(|i?|A|Qi\ii;|)(^-i)/^<c,Per(i?,Qi), 

and also the optimal constant for the Poincare inequality 

||/-(/)|lil*(Q,)<C/P/||Ll(Q,)' 

for / G W'HQi), I* = N/iN - 1) and (/) = /^^ /(x) dx. 

To state the main results we introduce some additional notation. Through- 
out the paper we use Qf and BUC{U) to denote respectively the space-time 
cube Q+ = Qi X [0, 1] C M^+i and the space of bounded uniformly continuous 
functions on U , and we write 

{V) = I V{x,t)dxdt and {xdivV) = / xdwV{x,t)dx dt . 
Jq+ Jq+ 

We have: 
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Theorem 1.1. Assume that V G C°'i(M^+i; M^) satisfies ([O]) and ([LS]) and 
that uq € C'^(M^) is bounded. There exists a positively homogeneous of degree 
one, Lipschitz continuous, convex Hamiltonian H : — > M such that, if S 
C''(M^ X [0, +oo)) and u € C^(M^ x [0, +cxd)) are the solutions to the initial value 
problems and respectively with initial datum uq, then, as e — t- 0, the 

u"^ 's converge locally uniformly in x [0, T] to u. Moreover, for all P £ M^, 

H{P) > \P\l^\l-ci\\dwVi;t)hN^Q^)) dt + {{V) + {xdivV),P). (1.6) 

Finally, the convex map P — H{P) — { {V) + {xdivV),P ) is coercive. 

For Lipschitz continuous initial datum uq, we can actually estimate the con- 
vergence rate as e — )• 0. We have: 

Theorem 1.2. Assume that uo G C°'i(M^) and let u'^,u £ BUC{W^ X [0,T]), 
for all T > 0, be respectively the solutions to and ()1.3p . Then, for all T > 0, 
there exists a positive constant C that depends only on T, N , V and the Lipschitz 
constant of uq, such that, for all {x,t) G x [0,T], 

\u{x,t) -u'{x,t)\ < Ce^/^ 

In the case that, for all x G M.^ and t G M, 

div^y(x,t) = 0, (1.7) 

we derive some additional properties of the function H. To simplify the statement 
we also assume, without any loss of generality (see Lemma l3. II below), that 

/ V{x,t)dx = ViGM. (1.8) 
JQi 

Then, according to Theorem II. H the averaged Hamiltonian H : — t- M 
satisfies, for all P G M^, 

H{P) > \P\. 

We establish here a necessary and sufficient condition to have the strict in- 
equality 

H{P) > \P\, 

in which case we have enhancement of the speed due to averaging. 

Recall that, since H is homogeneous of degree one, the level sets of u move 
with speed H{h) in the direction of the normal vector h = —Du/\Du\. Therefore, 
we refer to the situation H{P) > \P\ as "enhancement", because it implies that 
such velocity fields lead to faster propagation of interfaces compared to the case 
V = 0. 

First, we state the result in the case where V only depends on x. We have: 

Theorem 1.3. Assume that V G C°'^(M^) is -periodic, divF = and {V) = 

0, and let P G ]R^\{0}. Then H{P) = \P_\ if and only if, for all x G M^, 
(y{x),P) = 0. In particular, if N = 2, then H{P) = \P\ if and only if the stream 
function E associated to V is of the form E = E{{P,-)) for some : M — )■ 

1. e., V is a shear drift in the direction orthogonal to P. 

When V is also time dependent, the characterization of equality H{P) = \P\ 
is provided by 
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Theorem 1.4. Assume that, for allt G M, diva;y(-,t) = and Jq^ V{x,t)dx = 0, 
and fix P e M^\{0}. Then H{P) = \P\ if and only if there exists z G BVioc{V) 
suchthatz' > —\P\ in the sense of distribution and the function z{x,t) = + 
t) is Ij^^^ —periodic and satisfies, for all t € M, in the sense of distributions 

div((z(-,t) + (P,-))y(-,t)) =0 in R^. (1.9) 

We continue with some observations about these results. Theorem 11.31 yields 
that, if = 2, H{P) = \P\ and V is not constant, then P = {Pi,P2) must be a 
rational direction, i.e., either P2 = or P1/P2 G Q; since V is Z^— periodic. For 
Theorem 11.41 observe that, if z is not constant, then P/\P\ must be a rational 
vector. 

Also ()1.9p is equivalent to saying (see Lemma 13.21 below) that, for any fixed 
t > 0, the map x z{x,t) + {P, x) is constant along the flow of the differential 
equation X'{s) = V{X{s),t). 

We remark that it is possible to construct nontrivial examples of time-dependent 
flows for which H{P) = \P\. Indeed when N = 2 for any smooth, Z^-periodic 
{Ei,E2) such that Ei{0) = 0, let E{xi,X2,t) = Ei{xi + t)E2{x2), V = V^E 
and P = (1,0). Then H{P) = \P\ because the map z{s) = [s\ — s, where [s] 
stands for the integer part of s, satisfies the condition of Theorem 11.41 For more 
analysis and numerical computation of H for specific flow structures, we refer to 

[laiiaEoiEiiES]. 

The next result of the paper is about the long time behavior of the solution to 
(jl.ip with e = 1 and, in particular, the convergence, as t — >■ 00, of its zero level 
set to the Wulff-shape associated with the effective H, which is given by 

W = {y G : (P, y) + H{P) > for all P £ M^}. (1.10) 

We consider the initial value problem 

ut = \Du\ + {V {x,t) ,Du) in M^x(0,oo), 

(1.11) 

u = uo on X {0} , 
and set, for all t > 0, 

K{t) = {x eM.^ : u{x,t) > 0}. 
Recall that, in the language of front propagation (see, for example, [6J), the 
family of closed sets {K{t)t>o is solution of the front propagation problem 

v = l- {V{x,t),fi) 

starting from K{0) = Kq. 
We have: 

Theorem 1.5. Let Kq be a non-empty compact subset o/M^. There exist C > 
and T > such that, for all t >T, 

K{t) C {t + C)W. (1.12) 

Moreover, there exists a constant Cq > 0, independent of Kq, such that, if Kq 

contains a cube of side length Cq, then there exist C > and T > such that, 
for all t > T, 

{t-Ct^/^)W CK{t). (1.13) 
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We note that we do not know whether the size condition on Kq is actuahy 
necessary. 

The final result of the paper is about homogenization when V depends on x 
also at the integral scale, i.e., we are interested in the behavior as e — )• of the 
solutions to the initial value problem 



,e 



\Du'\ + {V{x,l,l),Du' 



m 



X (0,T) 



(1.14) 



M" = Uq on M"' X {0} , 

where V : x X M ^ is smooth, bounded, Z^-periodic with respect to 
the last two variables, i.e., for ah {x,y,s) G R^ x R^ x R, 



V{x, y + k,s + h) = V{x, y, s) , 
is divergence free in the fast variable, i.e., for all {x,y,s) € 

divyV{x,y,s) = , 
and satisfies, for all x € R^ the "smallness" condition 



R^ X R, 



V{x,y,s)dyds 



The homogenized initial value problem is 



Ut = H{x, Du) in 
u = Uq on R^. 



< 1 



X (0,T), 



(1.15) 



(1.16) 



(1.17) 



(1.18) 



We have: 



OfmN 



Theorem 1.6. Assume (fTTS]) . (fTTe]) and (fTTS]) . There exists H G C 
M^), which is positively homogeneous of degree one and convex with respect to the 
second variable, such that, for any initial condition uq G BUC{M.^), the solution 
to Iil.l4\ ) converges, as e ^ 0, locally uniformly in R^ X [0, T] to the solution 
u of (jl.lSp . Moreover H satisfies, for all {x,P) E 



H{x,P) > \P\ + 



V{x,y,s)dyds,P) 



(1.19) 



The paper is organized as follows. Theorem 1 1 . 1 1 and Theorem 11.21 are proved in 
Section[5J In SectionOwe prove Theorem 1 1 . 3 1 and Theorem ll.4l while Theoreir ll.51 
is proved in Section SI In Section O we prove an extension of Theorem 11.11 to the 
case where V = V{x,x/e,t/e) has large-scale spatial variation. The Appendix 
contains a proof of Lemma 12.31 which plays an important role in the proof of 
Theorem II. 1[ 

About the time this paper was completed, we learned about a similar but less 
general homogenization result obtained by different methods in [25] . In particular 
it is proved in [25] that homogenization takes place for time independent advection 
satisfying V = Vi + V2 with divVi = and IV2I < 1- 

Finally we remark that throughout the paper we will need some basic results 
from the theory of viscosity solutions, like comparison principles, representation 
formulae, etc.. All such results can be found, for instance, in [5]. 
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2. HOMOGENIZATION 

We begin with some preliminary discussion and results to set the necessary 
background for the proofs of Theorem 11.11 and Theorem II .21 First, we recall that 
for any A > and any P G M^, the "penalized" cell problem 

vx,t + Xvx = \Dvx + P\ + {V,Dvx + P) in (2.1) 

has a unique Z^"'"^-periodic solution vx E BUC{M.'^~^^), which is actually Holder 
continuous and satisfies, for all {x,t) G M^^"*^, the bound 

- X-^\P\{1 + llFlloo) < ^^A < A- V|(l + Halloo). (2.2) 

We also recall that, in the periodic setting, homogenization is equivalent to 
proving that the (AfA)'s converge uniformly in M^+\ as A — > 0, to some constant 
c(P) and that c G C°(M^). In this case, H{P) = c{P). 

In view ()2.2p . to prove the convergence of the Af^'s, we need to control their 
oscillations. We have: 

Lemma 2.1. For all P eM.^ and Xe (0, 1], 

osc(Ava) < C\P\X, (2.3) 

where 

C7 = 4(l + ||F|U)(2i/^a*~'iV + 3). 

Before we present the proof of Lemma [2. 11 which is one of the most important 
parts of the paper, we point out its main consequence in the next 

Corollary 2.2. Let C be the constant given by Lemma [Ol There exists some 
H{P) £ M such that 

\\Xvx-H{P)\\oo<C\P\X. 

Proof: The maps A — > A min vx and A — )• A max vx are respectively nonincreasing 
and nondecreasing. For the sake of completeness we present a formal proof, which 
can be easily made rigorous using viscosity solutions arguments. Since the two 
claims are proved similarly, we present details only for the former. 

To this end, for < A < /i, let (x, t) be a maximum of — vx- Then, at least 
formally, at {x,t), we have D^^ti^^ = Dx^t^x^ and 

Vt,,t + l^Vf, < \DVf, + P\ + {V, DVf, + P) 

and 

vx,t + Xvx > \Dvx + P\ + {V, Dvx + P) • 

It follows that 

fJ-Vf,{x,t) < Xvxix,t) 

Let {y, s) be a minimum point of Xvx- Then 

l^v^{y,s) < n{v^{x,t) - vx{x,t) +vxiy,s)) 

< Xvx{x,t) - fivxix,t) + nvx{y,s) 

< Xvx{x,t) + \{vxiy,s) - vxix,t)) < Xvx{y,s), 

and, hence, 

min fiv^ < min Xvx ■ 
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The above remark combined with Lemma |2 . 1 1 imphes that the Xv\s converge 
uniformly to some constant S{P) and that 

Aminf;^ < H(P) < Xma,xv\ . 

a 

We continue with the 
Proof of Lemma I2.lt Without any loss of generality, we may assume that 
V is smooth. Indeed, if the result holds for any smooth V, then it also holds by 
approximation for any V G C*^'^. 

Recalling (|2.2p . wx = Xvx satisfies, in the viscosity sense, 

wx,t-\Dwxix,t)\- {Vix,t),Dwxix,t)) >-CoX in R^+i , (2.4) 

where 

C7o = 2(l + ||y|U)|P| . 

It follows that {x,t) — > wx{x,t) + CqXI is a viscosity super-solution of the 
level-set initial value problem 

{i) zt = \Dz\ + {V,Dz) in R^x(0,oo), 

(2.5) 

(ii) z = wx on x {0} . 

The standard comparison of viscosity solutions then implies that, for all (x, t) € 

wx{x,t) + CoXt>z{x,t). (2.6) 

The next step is to understand the evolution of the perimeter of the level-sets 
of z. For this we need the following result, which is proved in the Appendix. 
We have: 

Lemma 2.3. Assume that V G Ci'i(M^+i) and let z G BUC{R^+'^) be a solu- 
tion of h2.5\) (i). Then, for any level 9 G ( inf z(-, 0), sup 2:(-, 0) ) such that 

{z(-, 0) = ^} = d{z{-, 0) > 0} = d{z{-, 0) < 9} and |{z(-, 0) = e}| = , 

and for all t > 0, we have 

{Z{;t) =9} = d{z{;t) >9} = d{z{-,t) < 9} 

as long as {z{-,t) < 0} 7^ 0. Moreover the sets {z{-,t) > 9} and {z{-,t) < 9} have 
locally of finite perimeter, and \{z{-,t) = 9}\ = 0. Finally, for any compactly 
supported tp G C^(M^), the maps t — > I{t) = J|^^ t)>e} ^{^)dx and t — > J{t) = 
I{z(- t)<e} '^i^)dx are absolutely continuous and satisfy, for almost all t > 0, 

'^I{t) = / ^{x){l-{V{x,t),v{x,t)))dn''-\x), 



dt 
and 



J{t) = - I ^{x){l + {V{x,t),u{x,t)))dn^-\x), 



dt J{z{-,t)=e} 
where z^(x,t) denotes in the former identity the measure theoretic outward unit 
normal to {z{-,t) > 9} at x ^ d{z[-,t) > 9}, while in the latter is the measure 
theoretic outward unit normal to {z{-,t) < 9} at x ^ d{z{-,t) < 0}. 
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Continuing with the ongoing proof, suppose that there exists G M with 
mfwxi-,0) <6 < supwxi-,0) , 

such that 

{wx{;O) = e} = d{wx{;O)>0} = d{wx{;O) <e} and \{wx{;0) = 9}] = , 
and such that 

\{x € Qi : wx{x,0) < e}\ < 1/2 . 

For ah t > 0, set 

p{t) = \{x e Qi : z{x,t) < e}\ . 

Let [0,T) be the maximal interval on which p{t) < 1/2 for any t € [0, T). Note 
that T > because p{0) < 1/2. We claim that, for ah <ti <t2 <T, 

p{t2) - Pih) < - r a(5)p(s)(^-i)/^ds . (2.7) 



Indeed fix a positive integer R and let Qn = (— , "f The space periodicity 
of z gives 

p{t) = R-^ \{x G Qr : z{x,t) < e}\ . 
For > small, let Xh e C°'i(M^; [0, 1]) be such that 

= 1 in Qr and Xh = in R^\QR+h , 
and, for any t € [0, T], set 



and note that. 



PR,h{t) = R ^ Xh{x)dx , 

J{z{-,t)<e} 



\Ym pR^h{t) = p{t) . 



It follows from Lemma 12.31 that, for almost all t G (0, T), 

^,PR,h{t) = -R-"" [ Xh{x){l + {V{x,t),u{x,t)))dn''-\x) . 

Moreover 

xh{x)dn''~\x) > n''-\{z{;t) =e}n Qr), 



j 

Ju 



i{z{;t)=e} 

and, in view of the spatial periodicity of z, 

n''-H{z{;t) = e}n Qr) > R''n''-\{z{;t) = d}n Qi) . 

The isoperimetric inequality in the box Qi and the fact that \{z{-,t) < 9} OQil < 
1/2 give 

n^'-HM; t) = e}n Qi) > - \{z{;t) <e}n Qii^^-^)/^ . 

ci 

Using once more the space periodicity of z we get 

\{z{.,t) <e}n Qi|(^-^)/^ ={R + i)-(^-i) \{z{.,t) <0}n Q«+i|(^-^)/^ 

>(i?(i? + i)-i)^-i(p«,,(t))(^-i)/^. 
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Combining all the above we obtain 



{z{;t)=e} CI 
Next we estimate the integral 



Xh{x){V{x,t),u{x,t))dn''-\x). 

i{zi;t)=e} 

For some constant k depending only on N we have: 

Xh{x){V{x,t),u{x,t))dn''~\x) 



{z{;t)=e} 



{z{;t)<e} 



div{xhV){x,t) dx 



< 



Xhix)divV{x,t)dx + ||-Dxh||oo||l^||oo \QR+h\QR\ 



i{z(-,t)<e} 
< ( / \dwV{x,t)f dx^/^i 



{z{-,t)<e} 



+ kR^-'\\V\\oo 



Hence, for almost all t £ (0,T), we get 



+ kR^-'\\V\\oo. 



d , , A { R 
j^PuAt) < -(- ( ^ 



N-l 



R+1 



k 



^ ||diVy(.,t)||^.(Q^))(p^„,(t))(^-l)/^ + -||F 

Integrating first over [^1,^2] and then letting /i — > and R — > +00 gives (|2.7p . 

Since a{t) > 0, p{t) is non increasing on [0,T). Hence T = +00. Integrating 
(|2.7p over (0, t) we obtain, for every t > 0, that 

r 1 f* 1^ 

/^(t) < [p^/^(0) - - a{s)ds^^ , 

where = max{s,0}. 

From the assumption p(0) < 1/2, it follows that 

and, hence, p = in [t* , 00). The continuity and the spatial periodicity of z then 
yield that 



z > 



m 



X [t*, 00). 



Let k be an integer in the interval [t*,t* + 1]. The space-time periodicity of 
wx and (|2.6p give 



inf inf wxixA) 

tG[0,l]a;€~" 



= inf inf w\{x,t) 

te[k,k+i\ xeK^ 

> inf inf zix,t) - CoX(t* + 2) 
te[k,k+i] x-eRJv 

> e-CoX{t* + 2). 



(2.8) 
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It follows that, if € M is such that 

\{x£Qi:wx{x,0)<e}\<l/2, (2.9) 

then 

inf inf wx(x,t) >e-CX (2.10) 

te[o,i] i^eM^ 

where 

C = Coit* + 2). 

Now we derive an upper bound. To this end suppose that € M with 
inf ti;;v(-, 0) < 9 < supw\{-,0) , 

such that 

{wx{; 0)=e} = d{wx{; 0)>e} = d{wx{; 0) < 9} and \{wx{-, 0) = 9}\ = , 
and such that 

\{xeQi:wx{x,0)>9}\<l/2. (2.11) 

The claim is that 

ma^wx{x,t) <e + CX, (2.12) 

x,t 

where C = 2\P\{l + \\V\\oo){2+N2-^/^ (a*)-^). Indeed, arguing by contradiction, 
we assume that 

maxwx{x, t) > 9 + CX . 

x,t 

Then, by continuity and periodicity of wx, there is some r £ [0, 1] such that 

\{x£Qi : wx{x,t) > 9 + CX}\> . 
Let z satisfy ()2.5l -(i')) with initial condition z = wx on M.^ x {r}. As before, we 



have, for all {x,t) G x [r, +oo), 

wx{x, t) > z{x, t) — CoX{t — r). 

Set 

p{t) = \{xeQi : z{x,t)>9 + CX}\ t>T, (2.13) 

and observe that p is continuous with p(r) > 0. Then, arguing as before, we find 
that, for all t < ti <t2, 

p{t2) - p{ti) > r «(^) 1 - P(5)})(^-'^/'^ ds . (2.14) 



Since a{t) > for all t > 0, it follows that p is nondecreasing on [r, +oo). We 
claim that there is some r<t* = r + l + iV/(2^/^a*) such that p > 1/2 for 
t > T. Indeed, otherwise, p < 1/2 on [r, t*] and integrating (|2.14p over [T,t*] 
gives 

Pit*) > (^ £ a{s) dsr > \ , 

which contradicts our assumption. Let now k be an integer in [T, T + 1]. The 
time-periodicity of wx yields 

1/2 < p{k) = \{x(^Qi : z{x,k)>9 + CX}\ 

< \{x G Qi : wx{x, k) + CoX{k - t) > 9 + CX}\ 
= \{xeQi : wx{x,0) + CoX{k-T)>9 + CX}\. 
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It follows from k < {t* + 1) that 

Co{k-T)<C and |{x G Qi : wx{x,0) > 9}\ > 1/2, 

which contradicts the definition of 6, and, hence, ()2.12p holds. 
Finally set 

^ = sup|6'GM: \{x eQi : wx{x,0) <e}\ . (2.15) 

In view of the above, using (j2.1(jp and (|2.12|) . we get 

mill min w\(xA) > 9 — CX (2.16) 

and 

max max w\(x, t) <6 + CX , (2-17) 

te[o,i] a;eM^ 

where C = 2|P|(1 + ||F||oo)(3 + A^2-i/^/a*). It follows that osc(tt;A) < 2CA, and, 
therefore, (ESI). 



□ 

We proceed with the 
Proof of Theorem [HU : Let H{P) be defined by Corollary O The fact that 
the map P — )• H{P) is positively homogeneous, convex and Lipschitz continuous 
follows easily from the properties of (|2.ip and the comparison principle of viscos- 
ity solutions. 

To prove (jl.6p . first we perturb (j2.ip by a vanishing viscosity, i.e., for r/ > 
we consider 

^l,t + K - = + P\ + Dvl + P) in (2.18) 

which has a unique Z^"^^— periodic solution G i?C/C(M^"*'^) which is at least 
in C"^(M^'''^) and converges uniformly, as ^ 0, to v\. 

Integrating (|2.18p over and using the periodicity, we find 

/ Xvldxdt= j \Dvl + P\dxdt+ [ {V,Dvl)dxdt + {{V),P) . 
JQt JQt JQt 



Set 



e(x) = (P,x), {vl{t))= vl{x,t)dx and (0=/ ^{x)dx. 

JQi JQi 

Since both V and v\ are Z^^-'^-periodic, for each t G [0, 1], we have 
{V{x,t),Dvl{x,t))dx 

Qi 

' {v1{x,t)-{vlmdiyV{x,t) dx 
Qi 

{vl{x,t) - {vl{t))+ax) - {0)dwV{x,t) dx+ 



Qi 



iC{x) - {0)divV{x,t)dx. 
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Applying, for each t G [0, 1], Holder's and Poincare's inequalities yields 

{vl{x,t) - {vl)t + ^{x) - {0)divV{x,t) dx 

Qi 

< if \vl{x,t)-{vl)t + ax)-m''^^''''^dx)^''-'^/''{[ |divy(x,t)|^dx 



< CI {[ \Dvl{x,t) + P\dx) \\divV{-,t) 
JQi 



\L^{Qi) ' 
and, since 

iiix) - {S))diYV{x,t)dxdt = { {xdivV) , P), 



Xvl{x,t)dxdt > 



Al-cj||divy(.,t)||^iV(Q^)) / \Dvl{x,t) + P\dxdt 
Jo JQi 



+ {{V) + {xdivV),P) . (2.19) 
Finally, in view of (|1.5p . we have, for all t > 0, 

l-Cl\\dwV{;t)\\LN^Q^)>0, 

while the periodicity of vx implies, again for all t > 0, that 



\P\ 



/ {Dvl{x,t) + P)dx < / \Dvl{x,t) + P\dx . 
JQi JQi 



Letting and then A -> in (|2.19p we obtain ()1.6p . 

□ 

Error estimates for the periodic homogenization of coercive Hamilton-Jacobi 
equations were obtained earlier in p[Qj. Although the proof of Theorem 11.21 is 
almost the same as the one of the analogous result in |10) . we present it here for 
the sake of completeness. To simplify the presentation below we denote by C 
constants that may change from line to line but depend only on uq, ||-Duo||oo and 
V. 

We have: 

Proof of Theorem [HI : For all (x, t, P) e x M x M^, set 

H{x,t,P) = \P\ + ly{x,t),P) . 

To avoid any technical difficulties due to the unboundedness of the domain, 
first we assume that uq is (Me)Z^— periodic for some positive integer M, which 
implies that and u are also (Me)Z^— periodic, and we obtain the estimate 
with constant independent of M. Then we use the finite speed of propagation 
property of the averaged initial value problem, to remove the restriction on uq. 

Let vx = vx{-, ■■,P) € BUC{R^+^) be the Z^+^-periodic solution to and 
recall that the map 

P — > Xvx{-, •; P) is (1 + ||y||oo)-Lipschitz continuous; (2.20) 

to simplify statements heretofore we say the / is L-Lipschitz continuous if it is 
Lipschitz continuous with constant at most L. 

Fix T > and consider $ : x [0, T] x x [0, T] ^ M given by 

^/ N F/ \ / N f X t x — y\ |x — (t — s)^ 
#(x, t, y,s) = u (x, t) - u{y, s) " '^^a ( -, -, -^p- j ^ ^ ds, 
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where /3 G (0, 1), A G (0, 1) and 6 > are to be chosen later. 

Since $ is periodic in the space variables, it has a maximum at some point 
{x,t,y,s). The main part of the proof consists in showing that either t = or 
1 = for a suitable choice of A and 5. 



We argue by contradiction and assume that t > and s > 0. Since u is 
Lipschitz continuous and ()2.20p holds, by standard arguments from the theory of 
viscosity solutions, we have 

1^ — y| - ( . e^^^^ 



< C 1 + 



A 



and 



\t - s\ 



< C . 



We claim that 



t - s 



< H 



* +C ^— + A . 



(2.21) 



(2.22) 



£^ J ' ~ \ X 

Indeed for q,/3 > small, let (xq, to-, r^, z^) be a maximum point of 
given by 

z — y\ |x — (t — s)^ 



^i{x,t,y,r,z) = u%x,t) - £vx[y,r, 



(ley - + \z- + \er - tp) - ^(|x - + (t - if). 



2eP 2e 

1 

-— ^, , ,,. ^, , ^ 2 

Since {x,t) is the unique maximum point of '^i{x,t,x/e,t/e,x), we have that 
{xa,ta, yaifa-, Za) Converges to (x, t, t/e, £) as a — >• 0, with 



1 



^0+ 2a 
while ()2.20p implies that 



lim — [\eya - Xa_\ + \Za - Xa\ + \£ra -ta\ ) = 



a A 
From the equation satisfied by u'^ we have 

+ ~ + /3(t, - i) 

e a 

^ TT ( Xa Xqi y Xa £y« , X^ -Zq , n/ 
< H ,— , -s \ \ h P{Xa -X, 

\ e e eP a a 

while from the equation satisfied by vx we also have 

ia , ^ / Zqc y \ ^ TT ( Xa ^Va , V 

\- Xvx {ya,ra, I > H ( yQ,r„, \ -g— 

Using the bound on the oscillation of the Ai^a's in Lemma |2.H we get 



(2.23) 
(2.24) 

(2.25) 



Za y \ , /^\ \Za y\ ^ j^r ( Xa ^Va , Z^ y\ ta STq, 

'+C;A 5 — >H[ya,ra, h 



Combining the above inequality with (|2.25|) and using the regularity of H gives 



^(^^) + CA^^ > ^ + /?(t„ - i) 



-C{ -3 \ h P\Xa -X\) 



-Ci 



Xa 

ya 

e 



+ 



a 



•^l\Xa ^ya\ ^ \Za 2/| 



a 
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Using (I2:23|) and (l2:2il) we now let a to get 
Recalling (|2.2ip we obtain 



and, since A G (0, 1), we finally get (|2.22|) . 
We now show that 



-i>H(^]-C^. (2.26) 



e 

To this end, for a, (3 > 0, we consider 

T / N -/ \ / X i X - z\ \x - yl'^ (i - s)'^ 

^2{y, s, z) = u{y, s) + evx I -, -, -j^ 1 + + — — — + Ss 

2a 2 ^'^ fi ^ I I ^' 
which has a minimum at some point {ya, Sa, Za)- Using (j2.20p and the fact 
that {y,s) is the unique minimum of the map (y,s) — )• ^2{y,s,y), we find that 
{ya, Sa, Za) couverges to (y, s, y) as a — )• 0, with 

\za-ya\^^e^ and lim = 0. 



a A a-i>o a 

Since li solves (11.311 we have 



(5 > H{ — /3[ya-y)) , 

e £P a 

and, in view of the Lipschitz continuity of H, letting a — )• 0^ gives ()2.26p . 

Combining (f2:22]l and ([2:26]) we obtain 

6 - C (^-^ + A^ < , 

which for A = ^(i^^)/^ and 5 = SCe'-"'^^^^/^ gives a contradiction. So either t = 
or s = 0. 

Next we estimate max^j^y^s ^(2;, 0, y, s) and max^j^t^j, <I>(x, t, y, 0). We have: 
max 0, y, s) 

x,y,s 

( E^~^ \x — yP 

< max I no(x) - uoCy) + Cs + C— — |x - y| —5 5s 

x,y,s [ A zeP 2e 

( \x — 2/ 1 ^ 

< max < C\y — x\ + Ce 



x,y 



2eP 



< C[e'^ + ej < Ce^, 
and, similarly. 



max$(x,t,y,0) < C e'. 

x,y,s 
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Therefore, for any {x,t) € x [0, T], we have 

u^x,t) -u{x,t) < 5t + Cel^ < CTe^^-^^l^ + Ce^ . 
Choosing /3 = 1/3 finally gives 

u\x,t) -u{x,t) < Ce^l^ . 
This reverse inequality can be obtained in the same way. 

□ 

3. ENHANCEMENT OF SPEED 

To prove Theorem 11.31 and Theorem 11.41 we need three results which we for- 
mulate as Lemmas next but present their proof at the end of the section. 
We begin with 



Lemma 3.1. Let V and H be as in Theorem li.il Then, for any Z— periodic 
c G C^(M;M^), the averaged Hamiltonian associated to V — c is H — {J^ c{s)ds, •). 

The second is 



Lemma 3.2. The divergence zero condition U.9\) is equivalent to the fact that, 
for any fixed time t, the map x — )■ z{x, t) + {P, x) is constant along the flow of 
the ode X'{s) = V{X{s),t). 

To state the final result we recall the notion of an e-mollifier. To this end, 
let (j) G Cf'(M^+^[0, 1]) be such that ^(0) = 1 and J^N+ifp = 1 and define 
the e-mollifier (pi; by (ps = e~^'^^^^(j){x/e,t/e). Then /jgjv+i (pe = 1 and, for any 
/ G LI^^{M.'^~^^), /e = / * is a smooth approximation of /. 

For the rest of the section we assume that, for each t G M, 

div^.y(-,t) = and / V{x,t)dx = 0; (3.1) 

JQi 

recall that the average zero condition is actually not a restriction, since we can 
always replace V by ^ ~ Jq-^ ^(^j t)dx a fact, which, in view of Lemma 3.1, simply 
adds a translation to the effective Hamiltonian. 
The final preliminary result is 

Lemma 3.3. Assume ()3.ip and fix P £ M^. There exits a bounded, Z^^^- 
periodic z G BVioJ^^ x M) such that, for all e > the smooth functions z^ = 
z * 's satisfy in M^+^ 

Ze,t + H{P) >\DZe + P\ 

+ [ {z{y,s) + {P,y)){DM--y,--s),V{y,s))dyds . ^^'^^ 

We proceed with Theorem 11.31 which is a straightforward consequence of The- 
orem [T31 We have: 

Proof of Theorem [113} If H{P) = \P\, then let 

{P,x) 

z{x,t) = z{—-y- — ^ t) and zi{x) = / z{x,s)ds, 
I "I Jo 

where z is given by Theorem 11.41 
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Since V is independent of x, zi also satisfies (|1.9p in the sense of distributions, 
and, since z is periodic zi is actually a constant. In this case (|1.9p reduced to 



{V{x),P) = for all x £R^. 

When N = 2 and {V) = there exists a Z^-periodic stream function E -."M? ^ 
M such that V = (-^, In this case, if P = (^1,^2) and q = (-Pi,P2), 

(y, P) = in m? becomes 



If q is an irrational direction, then the map t ^ x + tq\s dense in R^/Z^. Since 
E is constant along this trajectory, E is constant and therefore V is identically 
equal to 0. Otherwise, t ^ x + tq \s periodic and E has to be constant along 
this trajectory. This means that E = E{{-,P)) for some smooth periodic map 

: M — > M, and, hence, F is a shear advection. 

Conversely, if {V, P) = in M^, then vx = \P\/X is the unique solution to 
(j2.ip . and Xvx = |P| clearly converges uniformly to H[P) = |P|. 



We turn now to the 

Proof of Theorem [TH Assume that, for some P G M^\{0}, H{P) = P and 
let z be given by Lemma 13.31 We first prove that z is a function of only (P, x) 
and t. More precisely, we claim that there exits z E BVioci^) such that 



To this end, let (ps and z^ = z * (j)^ he as in Lemma l3.3t note that z^ is Z 
periodic. Then, for ah (x, t) G x M, 



{)={DE{x),q) . 



□ 



z = z{{P,x),t) and Zs{s,t) > —1 in the sense of distributions. 



(3.3) 




{z{y,s) + {P,y)){D(t)e{x - y,t - s),V{y, s))dyds dx = . 



(3.4) 



Indeed the periodicity of z and V gives 




iz{y,s) + {P,y)){D(l)e{x -y,t- s),V{y, s))dyds dx 




{z{y, s) + (P, y))^e{x -y,t- s){V{y, s), u^)dyds dn^-\x) 



(z{x - y,t - s) + {P,x - y))(l)e{y,s){V{x -y,t- s),Ux)d'H (x) dyds 



{P,x- y)(pe{y,s){V{x -y,t- s),v^)d'H (x) dyds. 



while (|3.ip yields 




(P,x - y)(t)e{y,s){V{x -y,t- s),v^)d'H (x) dyds 




and, hence, (|3.4p holds. 
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Next we integrate (|3.2p over Qi x (0,1). Using (j3.4p and the periodicity of 
Zsi-,t) we get 

"1 r rl 



dt = \P\ 



\P\=H{P)> f [ \Dz,{y,t) + P\dydt> [ [ {Dz,{y,t) + P)dy 
Jo JQi Jo JQi 

It follows that, for all {x,t) G there exists 6{x,t) > —1 such that 

Dzeix) = e{x,t)P. Thusze isoftheformz£(x,t) = Ze{{x,P),t), with : ^ R 
satisfying Z£^s(s,t) > —1. 

Passing to the limit e — )• 0, we also find that z = z{{x, P) ,t) for some map 
z G -BV;oc(M^,M) satisfying Zs{s,t) > —1 in the sense of distributions. Whence 
([3:3]) holds. 



Next we claim that z satisfies, in the sense of distributions, 

St{s,t) = Ss{s,t)\P\ inM^, (3.5) 
and that (|1.9p holds. Note that this proves the "if" part, since (|3.5p implies the 



existence of a map z G BVioc(^,^) such that 

z{s,t) = z{\P\~'^s + t) . 
Moreover we have z'{s) > —\P\ in the sense of distributions because Zs{s,t) > 

\p\ 



-1 in the same sense. Finally, z{x, t) = z i -^jpr + t \ is periodic in space and time. 



We continue with the proofs of ()3.5p and (jl.Op . If z is constant, then ()3.5p is 
obvious and ()1.9p just follows from ()3.2p when e — )• 0. 

Next we assume that z is not constant. In this case z = z{{-,P),t) is Z,^^^- 
periodic and not constant. Therefore P has to be a rational direction. So, up to 
a rational change of coordinates, we may assume without loss of generality that 
P = Oei for some > 0, while V is still Z^"*"^— periodic. 

Using the notation x = {xi,x') for each vector of with xi G M and x' G 
M^~-^, for a fixed (xi,i) G we integrate (|3.2p over the cube {xi} x Q'^ x {f}, 
where Q\ = (-1/2, l/2)^-\ and obtain 

Ze,tmx),t) + \P\ > (4,s((P,x),t) + l)|P| + 

iz{e{xi - yi),t- s) + 9{xi - yi)){D(j)eiy,s), V{x -y,t- s))dyds dx' . 

-1 

(3.6) 

It turns out that the last integral in the right-hand side of the above inequality 
vanishes. Indeed 



{D(t)e{y,s),V{x - y,t - s))dy'dx' = 

-1 

/ ^(y^^) f Vi{x -y,t- s)dx'dy' 
J^N^i dxi Jqi^ 

+ / (l)e{y,s)—^{x-y,t-s)dy'dx'. 

Jq{ Jr^-^ dxj 

The periodicity of V yields, for any j = 2, . . . , N , 

4'e{y,s) [ ^{x -y,t- s)dx'dy' = , 
Jq{ 9xj 
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while the divergence free condition and, again, the periodicity give 

- — / Vi{x-y,t-s)dx' = —^{x-y,t-s)dx' = -y^ -^{x-y,t-s)dx={). 
dxi Jq'^ Jq'^ dxr ^ Jq,^ dxj 

On the other hand, 

// Vi{x — y,t — s)dx' dxi = f Vi{x,t — s)dx = 0, 
-1/2 Jq\ Jqi 

and, hence, for ah xi G M, 

/ Vi{x-y,t-s)dx' = . 

Therefore 

s), V{x -y,t- s))dy'dx' = , 



which, going back to ()3.6p . proves that 

4,t((P,x),t) >5,,,((P,x),t)|P| in M^+^ 

Since Ze is periodic, integrating the above inequahty over (—1/2, 1/2) x (0, 1) 
shows that in fact it must be an equahty. Letting e ^ then gives (j3.5p . 

To prove ([Tg]), we first combine ([32]), and ([33]) to get, for all (x,t) G 

0> / {z{y,s) + {P,y)){D(t>e{x-y,t-s),V{y,s))dyds . 



Averaging over the cube Qi we see that, as a matter of fact, equality must 
hold for all (x,t) G M^+i. 

Integrating the resulting equality against any compactly supported smooth 
function -0 : — )• R we get 

ip{x){z{y,s) + {P,y)){Dcl)s{x -y,t- s),V{y,s))dyds dx = 0, 

Kiv+i 

and after integrating again by parts and letting e — )• 0, we obtain, for all i G R, 

{Di;{x),iz{x,t) + {P,x))Vix,t))dx = 0, 



which is exactly (jl.9p . 

To prove the "only if part let z be as claimed. Since z G -BV;oc(I^) I^)) we may 
assume without loss of generality that z, and, hence, z{x,t) = + t) are 

lower semi-continuous. 

We show next that z satisfies, in the viscosity sense, 

dtz + \P\>\Dz + P\ + {V,Dz + P) inR^+i. (3.7) 

To this end, let (j) he a, smooth test function such that z > (p with equality at 
{x,i). It follows from equality z{x,t) = z{^^py- + t) that D(j){x,t) = 6P/\P\ where 

Since z' > —\P\ in the sense of distribution, it follows that 9 > —\P\- Moreover 
recalling (jl.Op and Lemma 13.21 we have, for a fixed t, that the function x — )• 
V{x,t) + {P,x) is constant under the flow of the ode X'{s) = V{X{s),t). 
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Let now X be a solution with ^(0) = x and t = t. Then, for any s E M, 
(l){x,t) + {P,x) = z{x,t) + {P,x) = z{X{s),t) + {P,X{s)) > ^{X{s),i) + {P,X{s)) , 
and, therefore, 

= [^{X{s),t} + {X{s),P)] = {D(t>{x,i) + P,V{x,t)) . 

as Is^o 

Combining the above relations gives 

(l)t{x,i) + \P\ = \D(l){x,t) + P\ + {V{x,i),D^{x,i) + P) , 
i.e., z is a super-solution of ([37 



Then it is easy to check that z{x,t) = z{x,t) + + ||-z||oo is a super-solution 
of ()2.ip . It then follows from the comparison principle that 

vx < -Z + A-^l + ||z||oo . 

Recalling that Xvx converges uniformly, as A — ?■ 0, to H{P), we obtain H{P) < 
\P\. Since the reverse inequality always holds, the proof of the implication is 
complete. 

□ 

We continue with the proofs of the lemmas. We begin with the 
Proof of Lemma [STU Fix A > and P € , let vx be the solution of (pT]) 
and us recall that Xvx converges uniformly, as A — >■ 0, to H{P). 

Set c = Jq c{s)ds and consider wx € BUCiM.^^^) given by 

wx{x,t) = vxix - c{s)ds,t) - { {c{s) - c)ds + ^,P) + 2\\c\\oo\P\ ■ 
Jo Jo ^ 

It follows that Wx is a super-solution of 

zx,t + Xzx = \Dzx + P\ + {V -c,Dzx + P) mR^+\ (3.8) 

We only present a formal proof, which can be easily justified using viscosity 
solution arguments. To this end observe that it is immediate from the definition 
of Wx that 

wx,t + Au;a > -{Dvx, c{t)) + vx,t - (c - c, P) + Xvx - (c, P), 

while 

\Dwx + P\ + {V -c, Dwx + P) = \Dvx + P\ + {V -c, Dvx + P) ; 

the comparison principle now gives wx > zx where zx is the solution of (j3.8p . 

Since the A2;;^'s and the Xwxs converge uniformly, as A ^ 0, to the averaged 
Hamiltonian Hc{P) associated ioV — c and to S{P) — (c, P) respectively, we get 
Hc{P) < H{P) — {c, P). The opposite inequality is proved similarly by considering 
— c instead of c. 

□ 

We continue with the 

Proof of Lemma [SH : Let X^{-) be the solution of X'{s) = V{X{s),t) 
with initial condition x at time s = 0. Then, for any /i € M and Tp : C^{R^) the 
divergence zero property of V yields 

ilj{x){z{X,{h),t) + {P,X,{h)))dx = [ iP{X^{-h)){z{x,t) + {P,x))dx . 
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Therefore, in view of (jl.Op . 

-^^ I i^{x){z{X^{h),t) + {P,X^{h)))dx 

{Dilj{x),V{x,t)){z{x,t) + {P,x))dx = , 

Applying this last equality to the test function ip o Xx{—s) we get 

/ i;{x){z{X^{s + h),t) + {P,X,Xs + h)))dx 

^ il,{X^{-s)){z{X^{h),t) + {P,X^{h)))dx = Q. 



dh 



h-Q 



Hence / ip{x){z{X''^ {s),t) + {P, X^ {s)))dx is constant in time, which means 
that z{-,t) is constant along the flow. 

□ 

We conclude with the 
Proof of Lemma 13.31 : For A > 0, let vx be the solution of of (|2.ip and set 

zxix,t) = vx{x,t) - vx{0,0). 

It follows from Lemma [2. II and Corollary 12 . 21 that . for some C > independent 
of A, 

II^aIIoo < c\p\. 

Next we show that the zx^s are also bounded in BViod^^'^^)- Indeed, for 
a > 0, consider the Z^+^-periodic solution vx,a to 

vx,a,t + Xvx,a = aAvx,a + \Dvx,a + P\ + {V, Dvx,a + P) m , (3.9) 

which is at least in C^'^ and, moreover, converges uniformly, as a — ?> 0, to vx- 

Integrating (|3.9p over a cylinder of the form Qr x {—R,R) for some positive 
integer R, we obtain, using the periodicity, that 

/ Xvx,a> [ \Dvx,a + P\ 

JQrx(-R,R) JQiiX{-R,R) 

+ {V{x,t),Dvx,a{x,t) + P)dxdt . 

jQRxi-R,R) 



Since, in view of (j3.ip . 

/ {V{x,t),Dvx,a{x,t) + P)dx = 0, 



it follows that 



[ \Dvx,a + P\<2R^+^\\Xvx,a\\oo . (3.10) 

JQiix{-R,R) 
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Let ^> e C^iQR X {-R,R)-R^+^) be such that < 1 for all {x,t) £ 

iN+i^ From ([32]) and ([XTH]) we get 



/ vx,adwx,t'^dxdt = - [v\,a,t'^ + {Dvx,a, '^)]dxdt = 

JQrx{-R,R) JQiiX{~R,R) 

- [ [{^Vx,a - a^Vx,a " \Dvx,a + P\ - {V, Dvx,a + P))^ + {Dvx,a, ^)]dxdt 

JQrx{-R,R) 

< CR{P)\\Xvx,a\\oc + a vx,a^^dxdt, 

JQrx{-R,R) 

where Cr{P) depends only on N, R, ||l^||oo and P. 
Letting q — )• yields 

[ vxdwx,t^dxdt < Cr{P)\\Xvx\\oo < Cr{P){H{P) + C\) , (3.11) 

JQrx{-R,R) 

which in turn implies, in view of the assumptions on that the zx^s are bounded 
in L°° and in BVioc uniformly with respect to A. Hence the zx^s converge, up to 
a subsequence and in Lj^^, to some Z^'^^- periodic z G BViod^^^^)- 

Let e > 0, and as in statement of the lemma, set zx,£ = 4>e * zx, and 
z\,a,e = 4>e * zx,a and fix {x, t) E R^"^-*^. 

It follows from ([gTH]) that 

ZX,a,s,t + >^Zx,a,e + Awa,q(0, 0) > aAzx,a,e + \Dzx,a,e + P\ 

+ / 4>eix - y,t - s){V{y,s),Dzx,aiy,s) + P)dyds . 
Integrating by parts and using (jS.ip we find 

ZX,a,e,t + >^Zx,a,6 + At;A,a(0, 0) > aAzx,a,e + I^^A,q,£ + ^| 

+ / (zxAv^^) + iP^y))iP'^eix - y,t - s),V{y,s))dyds . 



Letting first a — and then A ^ gives (I3.2p . 

□ 



4. Convergence to the Wulff shape 

We begin recalling some important facts from the theory of front propagation. 
The first is (see, for instance, [6]), that the family of sets {K{t))t>o is independent 
of the choice of uq as long as Kq = {x € : uo{x) > 0}. The second (again 
see [B] and the references therein), which we will use repetitively in the sequel, 
is the following superposition principle of the geometric fiow. If {Kq)q^q is an 
arbitrary family of non-empty closed subsets of M^, with corresponding solution 

{K^ {t))t>o, then the solution starting from (JeeO given by (Ueee -^^(*))t>o- 
This can be seen either by using the control representation of the geometric flow 
or the stability and comparison properties of viscosity solutions. A consequence 
is the well known inclusion principle. If Kq C K'q are two non-empty closed 
subsets of M^, then the corresponding solutions {K{t))t>o and {K'{t))t>Q satisfy 
K{t) C K'{t) for ah t > 0. 

We also remark that the convergence, as i — > oo, of the level sets of geometric 
equations without spatio-temporal inhomogeneities was considered in [18j . The 
results of [18] do not, however, apply to the problem at hand. 
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The proof of Theorem 11.51 is long. We formulate two important steps as sepa- 
rate lemmas which we prove at the end of the section. 

The first is about some "controllability" estimates. We have: 

Lemma 4.1. Assume (|3.ip . There exist a positive integer hq and T > such 
that, for all x € M^, the solution {K{t))t>o of the front propagation problem 
starting from the set x + {k ^ : \k\ < uq} contains Qi{x) at time T. 

The second is about some growth property for fronts. 

Lemma 4.2. There exits a positive integer R and positive constants r and Ti 
such that, if the initial compact set Kq C contains a set of the form Qj^{k) 
for some k G Z^, then the solution {K{t))t>o of the front propagation problem 
starting from Kq satisfies, for all t >0, 

Qrt{k) <ZK{t + Ti) Vt>0. 

We continue with the 
Proof of Theorem 11.51 : We first show ()1.12p . It is well known (see [6] ) that 
the characteristic function '^K(t) of K{t) is a solution to the geometric equation 

dtu=\Du\ + {V,Du) in X (0, +oo) . (4.1) 

Fix next a direction v G with \v\ = 1. For A > 0, let v\ be the solution of 

vx,t + \vx = \Dvx + y\ + {V, Dvx + v) in M^+^ , 

and recall that there is some constant C independent of v and A such that 

osc{vx)<C and \\Xvx - H{v)\\^ < CX . 

It is immediate that 

z[x, t) = vx{x, t) - va(0) + {v, x) + (^(z^) + CX)t + C , 

with 

C = C — min {u, x) + 1 , 

y£Ko 

is a super-solution of (|4.ip . Since ()4.ip geometric, it is immediate that max(2;,0) 
is also a super-solution. Moreover, we clearly have max(2;(-, 0), 0) > Irq- 
The standard comparison gives, for ah {x,t) G X [0,-Foo), 

max(z(x,t),0) > lK{t){x) , 

which, in turn, implies that, for all t >0, 

K{t) C {x G : z{x,t) > 1} . 

Recalling that H{jy) > = 1, we find, for all {x,t) G x [l,oo), 

z{x,t) < {iy,x) + {H{iy) + CX)t + C + C < {u,x) + H{iy){{l + CX)t + C + C) , 

and, therefore, for all t > 0, 

K{t) c{x£R^ : {u, x) + H{v){{l + CX)t + C7 + C - 1) > 0} Vt > 1 . 

Letting A — ?• 0, we get, for all t > and a new positive constant C, 

K{t) C {x G : {v, x) + H{v){t + C) > 0} > 1 . 

Taking the intersection of the right-hand side over all v we obtain, by the 
definition of W, that, for all t > 0, 

K{t) C (t + C)W . 



(4.3) 
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The proof of (jl.lSp is more intricate. To this end, let R, r and Ti be defined 
by Lemma [^2] and let [K{t))t>o be the solution of the front propagation problem 
starting at Kq C M^, a compact set which contains Qj^{k) for some k > Z^. 
Then, for all t > 0, 

BrtCKit + Ti). (4.2) 
Recall that u^{x,t) = lx(t/e)ix/^) is the solution to 

uf = \Du''\ + {V{x/e,t/e),Du'') in M^x(0,oo), 
n^(x,0) = lKo{x/e). 

and, in view of (gj]), for ah {x,t) G 

u'{x,t) > 1b^(^/^_^^j(x/£) = 1b^(^_^^^)(x) . 

Fix next 5 € (0, 1) such that Ti + 5 /{re) is an integer and, for ws{x) = {5 — 
|x|) VO, let be the solution of (|4.3p with initial datum w^- Then we know from 
Theorem 11.21 that there exists a constant C > such that, for all t G (0, 1), 

where is the solution of the homogenized problem 

ws^t = H{Dws) in x (0, 1) , 
ws{-,0) = ws on . 

Note that the constant C is independent of 5 because the ws^s have Lipschitz 
constant which are bounded uniformly in 6. 

From the Lax-Oleinik formula, iBs is given, for (x, t) € x [0, oo) by 

wsixjt) = sup ws{x — ty) . 

y&'W 

Since ws < Isg < u^{-,eTi + d/r), it follows, from the time-periodicity of V 
and the choice of 6, that, for all (x,t) € x [0, 1]), 

wsix,t) <u''{x,t + eTi + 5/r) . 

Hence, for all t G [0, 1], 

{wl{-,t) > 6/2} c {u'{-,t + eTi + 5/r) > 6/2}, 

while 

M(-, i) > 5/2} D {wsi-,t) > 5/2 + Ce^/^} D { sup wsi- - ty) > 5/2 + Ce^^^} . 

yew 

If we choose 5/2 - Ce^/'^ > 0, then 

{ sup wsi- - ty) > 5/2 + Ce^/^} D {x : sup {5 - \x - ty\) > 6/2 + Ce^/^} D tW . 
yew yew 

Therefore, for all t E [0,1], 

eK{{t + eTi+6/r)/e) = {u%-,t + eTi + S/r) > 6/2} D tW , 
i.e., for ah t G [0, 1/e], 

K{t + Ti + 6/{re)) D tW . (4.4) 

Finally, for t sufficiently large, choose e = \/{t — ACt^^^ /r) and 6 = {n — Ti)re 
where n is the integer part of ACt^^^ /r + 1. Then n = Ti + 6 /{re) is an integer. 
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5/2 — Ce^/^ is positive and, applying inclusion (j4.4p to t — n which belongs to 
[0, we get 

K{t) D {t- Ct'^/^)W 

for some new constant C. 

□ 

We conclude the section with the proofs of the two lemmas used in the proof. 
We have: 

Proof of Lemma 14.11 : Fix x G and let {K{t))t>o be the solution of 
the front propagation Vx^t = 1 — y{x,t) starting from Kq = {x} + Z^, i.e., 
K{t) = {y € : u{y, t) = 0} where u is the solution to 

dtu = \Du\ + {V, Du) in x (0, oo) 

u = —cLkq on R^, 

where dxo is the distance function to the set Kq. 

Then, for each t > 0, the set K{t) is Z^-periodic and a non empty interior 
(because it has an interior ball property, as recalled in Appendix). So p{t) = 
\K{t) r\Qi{x)\ is positive for positive time and, following the computation in the 
proof of Theorem 11.11 it satisfies for all t2> ti> 0, 

p{t2) - p{ti) > - r (min{p(i), 1 - pit)})^''-^^/'' dt . 

CI Jti 

Hence there exists a time T depending only on N such that \K(T)r]Qi{x)\ = 1. 
This means that Qi{x) C K{T). 

It follows from the finite speed of propagation, that there exits a positive integer 
no such that the solution K{t) starting from {x + k £ : |A;| < no} coincides 
with K{t) on Qi{x) x [0,r]. Then Qi{x) C K{T). 

□ 

Proof of Lemma 14.21 : Consider the solution u to 
ut = H{Du) in M^x(0,oo), 
t((x,0) = V (-1) on X {0} . 

Since, H{p) > \p\, given < 9i < Oq < 1 and 6 > small, there exists t G (0, 1) 
such that 

M-,0) > -9o} + B{0, 6) C {u{-,t) > -9i} . 
The fact that the solution of 

dtu' = \Du'\ + {V{^,^J,Du') in M^x(0,oo), 
n^(3;,0) = (-|x|) V (-1) on x {0} , 

converges, as e — )• 0, locally uniformly to u, yields an e € (0, 5/4) such that 

K(-,0) > -eo} + B{0,^-) c K(-,t) > -9o} . (4.5) 

Next fix no and T as in Lemma I4.H choose e such that Tq = - is an integer 
and 

Q„o+i(0)ci?(0,^), (4.6) 

and set 

Ko = {x : u'{ex,0) > -9o} . 
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The solution of the front propagation problem starting from Kq is given by 
K{t) = {x : u'^{ex,et) > —9o}. From ()4.5p we have 

while from (|4.6|) . for any k € with \k\ < no, we have Kq + k C K(Tq). 

The periodicity of V also implies that the solution of the front propagation 
problem starting from Kq + k is just K{t) + k while the solution starting from 
K{To) is Kit + To). 

From the inclusion principle we get, for all k E with jA;| < no and all t > 0, 

Kit) + kcKit + To) . (4.7) 

Then Lemma l4 . 1 1 implies that, for all t >0, 

Kit) + Qi C Kit + To + T) . 

In particular, by induction, we get, for all positive integers n and all t >0, 

Kit) + QnC kit + n(To + T)) . 

Choose a positive integer M such that, for all t S [0, To + T], Kit) C Qm- 
Then, for all positive integers n such that n > M and all t € [0, To + T], 

Qn-M C Kit) + QnC kit + n(ro + T)) , 
and, hence, there exist r > and Ti > such that, for all t >0, 

Qrt C kit + Ti) . 

Finally choose a positive integer R such that Kq C Qh- Then, for any compact 
initial set Kq such that Qj^ik) C Ko for some k S Z^, we have ko + k C Kq. 
Therefore the solution of the front propagation problem Kit) starting from Ko 
satisfies, for all t > 0, 

Qrtik) C kit + Ti) + kC Kit + Ti) . 

□ 

5. HOMOGENIZATION FOR X-DEPENDENT VELOCITIES AT SCALE ONE 

Before we begin the proof, we remark that, since we are only able to prove 
that H is continuous with respect to the x-variable, uniqueness of the solution 
to equation (jl.l8p could be an issue. This is not, however, the case because, in 
view of (|1.17p and (|1.19p . H is coercive. Note that this is the reason we do not 
consider V that also depends on a slow time variable, because then the coercivity 
of the averaged Hamiltonian would no longer ensure a comparison principle for 
the limit problem. 

We continue with the 

Proof of TheoremlUBt For any fixed (x, P) £R^x M^, let ^ = v^'^'iy, t) 
be the solution to 

+ Auf ^' = \Dv^''' + P\ + {Vix, y, t), Dv^'"" + P) in x M . 
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r'"-^r"-iioo<CMiPi '"^"r^' (5.2) 



From Lemma l2.1l and Corollary 12 . 2 1 we know that there exist H{x,P) and, for 
any M > 0, a constant Cm > independent of A and P, such that, for all x € M.^ 
such that |x| < M, 

\\H{x,P) - Xv^'"" Woo <C\P\X . (5.1) 

Arguing as in the proof of Theorem 11.11 one easily checks that H is positively 
homogeneous of degree one and convex in P, and that ()1.19p holds. 

To complete the proof, it only remains to show that H is continuous in (x, P). 
To this end, observe that the standard comparison arguments imply H is {1 + 
||F||oo)-Lipschitz continuous in P. 

Fix next M > 0, xi,X2 € with |xi|, |x2| < M. Once again the standard 
comparison of viscosity solutions (see [5] ) gives 

a" 

where L = sup|^.|<M \\Dy^sV{x, •, •)||oo and Cm depends only on M. 

Combining ()5.ip with ()5.2p . we obtain, for all xi, X2 € such that |xi|, |x2| < 
M and |xi — X2I < 1, 

|^(X1,P)-^(X2,P)| < Cm\PH\xi - X2\) 

with (the modulus) co given, for r £ (0,0], by 

Lo(r) = inf (A- V^/^ + A) . 
Ae(o,i] 

The proof of the continuity of H is now complete. 

□ 

6. Appendix 

We present here the 

Proof of Lemma [23t Fix (9 e ( inf z(-, 0), sup z(-, 0) ) and let Kg{t) = {z{-,t) > 
9}. Since z is a solution to ()2.5p -(i). Kglt) is given by 

_ { ^ . '■ [0, — >■ absolutely continuous such that ^{t) = x. 



• z{m,0)>0, |e'(s) + y(e(s),s)| <la.e. sG(0,t) 
i.e., Ko^t) is the reachable set for the controlled system 

where /(x, s,a) = a + V{x, s) and a G are such that \a\ < 1. 

Since V is of class C^'^, it follows from [9J that, for any < r < T, there exits 
a constant r = r{T,T) > such that, for all t € [t, T], the set Kg{t) has the 
interior ball property of radius r, i.e., 

for all X € dKg{t) there exists y ^ such that x € B{y,r) C Kg{t), 

where B{y,r) stands for the closed ball of radius r centered at y, and, hence, for 
all t > 0, KQ{t) is a set of finite perimeter (see, for instance, [2]). 
Note that [9] deals with time independent dynamics, but the proofs can be easily 
adapted to the time-dependent framework considered here. 

Next set d{x,t) = dKg{t){x)- It follows that d is Lipschitz continuous in (x,t) 
and, moreover, 

dt = -l + {V{x - dDd, t) , Dd{x, t)) a.e in {d > 0} . (6.1) 
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It is, of course, clear that x d{x, t) is 1— Lipschitz continuous. The Ijy |loo + l- 
Lipschitz continuity of the map t — )■ d{x, t) comes from the above representation 
formula of Kg^f. 

To check (j6.ip . recall that, since the map P — >■ |-P| + {V^P) is convex, z is 
a sub-solution of ()2.5p (i) with an equality (see [8] and [E]), i.e., for any test 
function ip and any local maximum point (x, t) of z — ip, we have 

ipt = \Dip\ + {V,D^) . 

The invariance property of the geometric equation (|2.5p -(i) (see [6]) then im- 
plies that the map (x, t) — lKg(t){x) is also a sub-solution of ()2.5p (i) with equality. 

Assume now that d is differentiable at some point (x, i) and let y be the unique 
projection of x onto Kg{t) and note that w{y, t) := '^Kg(t) iu) + d{y + x — y,t) has 
a local maximum at {y,t). Indeed, if y ^ Kg{t), then, for all {y,t) sufficiently 
close to {y,t), 

w{y, t) = d{y + x-y,t) <l + d{x, t) = w{y,t), 

and, if y G Kg^t, then 

w{y, t) = 1 + d{y + X - y,t) < I + \y + X - y - y\ = 1 + \x - y\ = w{y, t) . 

Since {x,t) lxg(t){x) is a sub-solution of (|2.5p (i) with equality and |Z)d(x, t)| = 
1 we obtain 

-dt{x,t) = 1 + {V{y,i) , -Dd{x,t)) , 
and ()6.ip holds because y = x — d{x, t}Dd{x,t}. 

Next we assume that 6 is such that d{z{-,0) > 0} = {z{-,0) = 9}, which is 
true for almost all 6 £ (inf z{-, 0), sup z{-, 0)). It then follows from [6j that, for all 
t > 0, d{z{-,t) >9} = {z{-,t) = 9} for any t > 0. 

We now prove that, for ah (p G C^{R'^) and all t G (0,r), 

/ ip{x)dx— (p{x)dx = / / ip{x){l—{V{x,s),i'{x,s)))dT-L^^^{x)ds. 

JKg(t) JKeiO) Jo JdKg{s) 

(6.2) 

To this end, for h > small, let (^^ iM. —?■ [0, 1] be such that Ch{p) = 1 if P ^ 0, 
Chip) = 1 - if p G [0, h] and Chip) = if p > /i. 

Multiplying (j6.ip by ifC'^id) (which makes sense because the sets {d(-,s) = 
0} = dKgit) and {(i(-,s) = h\ have a zero measure since Kgit) and {(i(-,s) < h} 
have the interior ball property) and integrating over x (0, t) gives 

ipix)Chidix,t))dx - / ipix)Chidix,Q))dx 
Jm 

* -1 



. . (pix)i—l + {Vix — dix, s)Ddix, s), s), Ddix, s))) dx ds . 

/o J{0<d{-,s)<h} 

(6.3) 

Note that 



lim [ / (p{x)Ch{dix,t))dx— I (p{x)Chidix,Q))dx] = I ipix)dx— ipix)dx 

Jm Jm JKg{t) JKg{0) 
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We now concentrate on the right hand side of (j6.3p . Since \Dd\ = 1 a.e. in 
{d > 0}, the co-area formula impHes 



h 



ip{x){—l + {V{x — d{x, s)Dd{x, s), s),Dd{x, s))) dx ds 

{0<d{-,s)<h} 



t ^ fh 



h 



^{x){l - {V{x - av"{x, s), s), v"{x, s))) dn^-^{x) da, ds 

J{d{-,s)=a} 



where v"{x, s) is the measure theoretic outward unit normal to the set s) < 
o"}, which has finite perimeter since it satisfies the interior ball property with 
radius r + a. 

In order to complete the proof of ()6.2p we just need to use the following 

Lemma 6.1. Let E he a closed subset ofM.^ with the interior hall property of 
radius r > 0. Then, for all compactly supported in x ^ G C(M^ x S^~^), 

lim / <l>{x,i^''{x))dn^-^{x) = [ <^{x,i^{x))dn^-^{x) , 

'^^^J{dEi-)=(T} JdE 

where dsix) stands for the distance of x to E and i'^{x) (resp. i^{x)) is the 
measure theoretic outward unit normal to {dE{-) < cr} (resp. to E) at x £ dE. 

Proof: Set Eo- = {dE{-) < cr} and denote by Ha the projection of dE^ onto 
dE. It is known that IIo- is uniquely defined for 'H^~^ —a,.e. point x G dE^. Let 
= n^-^ [dEa and jla = UJi^a, /io = /io = n^'^ [dE. 

The first claim is that fia is absolutely continuous with respect to p,Q. For this 
let us first recall that, since E has the interior ball property of radius r, the map 
n~"^ is well-defined on := {x G dE ; 3y £ dE^ with |x — y| = a} and that 
n~-^ is Lipschitz continuous, with constant at most {r + a)/r, on F^ (see [2j). So, 
if Z is a Borel subset of dE, then 

fi^(z) = ^^-i(n;i(z)) < Lip(n^i)?^^^i(z) = ^/2o(^) . 

In particular, fi^ is absolutely continuous with respect to fio and, if f^ = 

dfiQ 

then /o- is bounded by 7^^^^— a.e. in dE. 
Next we note that, for any ^> € C(M^ x 5^"^), 

/ ^x,u%x))dn^'\x) = f ^y + au{y),v{y))fa{y)d'H''-'{y) . (6.4) 
JdE^ JdE 

Indeed, the definitions of /Xg- and /o-, give 

<^{y + av{y),u{y))fa{y)d'H''-\y) = [ ^y + a,,iy),u{y))d{UM{y) 
dE JdE 

$(n^(x) + au{Ua{x)), zv(n^(x)))(i/x^(x), 

dE^ 

which implies (|6.4|) because ^{Ilaix)) = y"{x) and li^ix) + cyv{Iia{x)) = x. 

In view of (|6.4|) . to complete the proof, it only remains to show that the {faYs 
converge, as cr —> 0, to 1 in L^{dE,'H^^^). 
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Applying (j6.4p with $ = 1 gives 
Per(S,) = 



dE 



Since < in dE, the lower-continuity of the perimeter implies 

that limo-^o Per(£^CT) = Per(£'). Using again the inequality < r~^r + a which 
holds 'H^~^—a.e. in dE, we obtain, in the limit a ^ 0, 



\l-Uy)\dn''-\y) < 



dE 



1 



r + a 



Per (OS) + 



dE 



N-l 



(y) ^ 0. 

□ 
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